Observing the fluctuating stripes in high T c superconductors 
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Resting on the gauge theory of topological quantum melting in 2+1 dimensions, we predict that 
a superconductor characterized by crystalline correlations on a length scale large compared to the 
lattice carries a new collective mode: the massive shear photon. This mode is visible in the electro- 
dynamic response and the ramification of our theory is that electron energy loss spectroscopy can 
be employed to prove or disprove the existence of dynamical stripes in cuprate superconductors. 



The explanation of conventional superconductivity by 
the BCS theory belongs to the highlights of twentieth 
century physics. However, it appears that BCS is too 
primitive to address the high T c superconductivity found 
in copper oxides. It is well established Q, 0, 0] that the 
superconducting liquid found in the cuprates is in com- 
petition with the so-called stripe phase |^ , which can be 
viewed as a special kind of electron crystal. By studying 
the behavior of the spin P, 0, and phonon Q sys- 
tem with inelastic neutron scattering an intriguing case 
emerged that these stripes survive in the form of tran- 
sient correlations in the liquid up to rather large length 
(~ 10 nm) and time (~ 10 meV) scales. This affair is 
still quite controversial because the data can also be ex- 
plained in other ways 0, 0, [HJ and there is a need for a 
unique, experimentally accessible signature for the 'quan- 
tum stripes' ^3- -"- n this Letter we present such a feature 
which we identified by studying the quantum field theory 
describing the topological quantum melting of crystals 
[lil llii . ITS . Hij | . We claim that there should be an extra 
propagating mode in the 'nearly ordered' superconduc- 
tor, reflecting directly the presence of the transient stripe 
correlations. 

One can already anticipate the presence of the extra 
mode by a simple quantum hydrodynamical considera- 
tion. Shear rigidity, the ability of the medium to re- 
spond reactively on a force acting in opposite directions 
on opposite sides of the medium, is uniquely associated 
with translational symmetry breaking, i.e., crystalline or- 
der. In a classical fluid one finds a viscous response on 
shear, but in the dissipationless superfluids there is no 
response to shear at all: these carry only compressional 
rigidity, turning into magnetic (Meissner) and electri- 
cal screening in the superconductor [16j. However, the 
'nearly ordered' superconductor should still have the ca- 
pacity to mediate shear forces at short distances and in 
analogy with the London penetration depth associated 
with magnetic forces, we call the characteristic length 
over which shear sources can make their presence felt 
in the medium the 'shear penetration depth' A5. The 
natural velocity scale of reactive shear corresponds to 
the transversal phonon velocity of the crystal ct, and 



assuming a coherent quantum dynamics, one anticipates 
the presence of 'photon-like' excitation in the system car- 
rying the shear stress, characterized by an energy gap at 
zero momentum, 57 = ct/\s : the 'massive shear pho- 
ton'. We claim that such a mode has to exist in high 
Tc superconductors when the dynamical stripe interpre- 
tation is correct, which implies that As 3> a, the lattice 
constant. Moreover, a controlled mathematical descrip- 
tion is available in this limit 0, ll^l IbH ITf| , based on 
modern developments in quantum field theory [l?! Il8l | , 
making it possible to arrive at a detailed, quantitative 
prediction for an observable quantity: the massive shear 
photon should give rise to a pole in the electron-energy 
loss function at small momenta and energies (Fig. 1), re- 
vealing its identity through the characteristic dependence 
of its pole strength on momentum (Fig. 1, inset). 

The applicability of the field theory requires some 
other conditions to be satisfied which are in fact im- 
plied by the assumption of a large crystalline correlation 
length. This means that the system has to be close to a 
continuous quantum phase transition between the crys- 
tal and the superconductor, and this implies in turn: (a) 
The liquid and the solid have to be microscopically sim- 
ilar and since the superconductor is a bosonic entity the 
crystal is also formed from bosons. It is about preformed 
pairs which can either form stripes or the superconductor 
|l9| . (b) The transition from the crystal to the supercon- 
ductor is first order, and to make possible a large shear 
length it appears necessary to assume that the supercon- 
ductor is at the same time a quantum liquid crystal of 
either the smectic or nematic kind |20j . It turns out that 
with regard to the shear stress photon it does not make 
much of a difference and we will present here results 
for the simpler nematic superconductor. 

The theory can be viewed as the quantum general- 
ization of the famous Nelson-Halperin- Young theory 
of classical melting in two dimensions, based on the no- 
tion that the liquid crystal can be viewed as a crystal 
where topological defects (dislocations) have proliferated. 
The quantum version rests on a Kramers- Wannier dual- 
ity arising in elasticity theory, first explored by Kleinert 
in the 1980's in the context of the classical problem in 
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three dimensions 1171. and only recently implemented on 
the quantum level |l2( . This earlier work contained some 
technical flaws (ignoring the relativistic character of the 
dislocation condensate, a faulty gauge fix) obscuring the 
view on the electrodynamic response. We present here a 
summary of the main steps of the correct derivation, and 
we refer for further details to Ref.'s 0HE1. 

In the 'orderly' limit As/a — > oo, the constituent 
bosons of the gaseous limit are no longer relevant and 
instead the work is done by the collective degrees of free- 
dom: the phonons and the topological defects. The lat- 
ter correspond to dislocations and disclinations, restoring 
the translational and rotational symmetry, respectively. 
In this language, liquid crystals are states where dislo- 
cations have proliferated, restoring translations, while 
disclinations are still massive excitations such that rota- 
tional symmetry remains broken [l5L lljj . The phonons of 
a ('Wigner') crystal of charged bosons in 2+1 Euclidean 
space-time dimensions are described by the following La- 
grangian density, following |l2l | from the gradient expan- 
sion in terms of the displacement fields u a , 

J? = ^u a C^ ab d,u b + iA a ^u a + \F^F^, (1) 

where C^ ua b is a short hand for the tensor of elastic 
moduli, now including the kinetic energy (p/2)(d T u a ) 2 
(Greek indices label space-imaginary time directions ~ 
x,y,T, Latin indices space directions ~ x, y; notice that 
C TTa b = p5 ab where p is the mass density). For simplicity 
we limit ourselves to isotropic elasticity characterized by 
just a shear and compression modulus, p and k respec- 
tively, which are related by the Poisson ratio k — p(l + 
f)/(l — v). This defines the velocities ct = \fp~Jp, Cl = 
c TV /2/(l - v), c K = c T a/(1 + f)/(l - v), corresponding 
to the transversal and longitudinal phonon velocities of 
the crystal, while c K is the compressional (real sound) ve- 
locity. The last term is the usual Maxwell term describing 
the electromagnetic fields while the electromagnetic vec- 
tor potentials (A^) couple to the elastic strains via the 
effective combination A^ = (n e e*) L4 r 0p a — A a S TtJ ] (e* 
is the microscopic electrical charge, n e is the density of 
charged particles). The scales associated with the elec- 



tromagnetism are the plasma frequency u> p = y/4r< 
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and the Debye electrical screening length A e = ct/ uj p and 
Debye momentum q e = 1/A e . 

The crucial insight, due to Kleinert, is to turn this 
into a 'stress gauge theory' (l7|. The first step is to use 
strain-stress duality, such that Eq. JIJ becomes, 
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in terms of the stress fields <r", with the er^'s hav- 
ing the status of canonical momenta. Notice that the 
antisymmetric components of the spatial stress vanish, 
a? — crjf = (Ehrenfest constraint). The key is that this 



is augmented by the law expressing the conservation of 
stress dfj,a^ — and this can be imposed in 2+1D by 
expressing the stresses in terms of the stress gauge fields, 



(3) 



The magic of these stress-gauge fields is, that like in vor- 
tex duality [2(|, the sources of these stress gauge fields 
correspond to the non-integrable displacement field con- 
figurations |l7|. 
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corresponding to the dislocation currents which can be 
factorized as J° = b a world-lines of dislocations with 
Burgers vector b = (b x ,b y ). This implies in turn that 
B'tJ® — > (b a B^) Hence, we arrive at a description of 
elasticity which is a close sibling of electromagnetism: the 
stress gauge fields ('stress photons') express the capacity 
of the elastic medium to propagate elastic forces, and the 
dislocations take the role of charged sources. 

In the 2+1 space-time dimensions of relevance to 
cuprates, the miracle is that the quantum-nematic crys- 
tal can be equally well viewed as a 'dual supercon- 
ductor' 0, ^3 1 m close analogy with the superfluid- 
superconductor duality in 2+1 dimensional phase dy- 
namics [HI, 13 113 ■ In the path-integral language, the 
'dual nematic shear superconductor' corresponds to a 
tangle of dislocation worldlines in space-time. As we ex- 
plained elsewhere [T^i these describe either smectic or 
nematic quantum liquid crystals depending on the orien- 
tational ordering of the Burgers vectors: in the smectic, 
dislocations condense with their Burgers vectors oriented 
in one particular direction while in the nematic they con- 
dense with equal probability for their Burgers vectors to 
be oriented along all allowed directions [Tg. As we are 
dealing with isotropic elasticity, we focus on the quan- 
tum hexatic, derived from the hexagonal crystal where 
the Burgers vectors are equally distributed along the six 
directions associated with the rotational Dq symmetry. 

The dislocation condensate is gauged by the stress 
gauge fields, and Eq. Q implies a covariant derivative 
structure of the form |(c^ — ib a B c ^)^\ 2 where \1/ is the or- 
der parameter of the dislocation condensate; as we over- 
looked in the earlier work 0] , it is vital to keep here the 
time derivatives given the relativistic nature of the dislo- 
cation condensate [2(|. Elsewhere, we will present in de- 
tail the averaging procedure leading to the correct Higgs 
term in the stress-sector 0,E3- We find a 'bare' Higgs 
term for the quantum hexatic J^H.bare = j\^q\ 2 B^B^, 
where ^0 is the exp ectation value of the disorder field. 
Both the glide |25j and Ehrenfest constraints have to 
be imposed on the bare Higgs term. The glide con- 
straint is implemented via a Lagrange multiplier term 
-2 "glide — iA(JJf — Jy ). This operation removes the com- 
pressional stress a% +qff from the Higgs term, having also 
the effect of turning [l2|j LLJ> the second sound velocity 
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of the condensate into the 'glide velocity' c g — ct/V%- 
Subsequently, the Ehrenfest constrained is implemented 
by eliminating the antisymmetric stresses crjf — a* from 
the Higgs term. It is particularly convenient to represent 
the resulting Higgs term in a gauge invariant way |l3| . 
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where we defined the 'shear Higgs mass' as fi = IvJ/ol^/M- 
This result is revealing: by counting derivatives it is 
immediately obvious that both the longitudinal (o% — <Jy) 
and transversal (a^, + a* ) shear stresses acquire a Higgs 
mass il, while compressional stress is left unaffected. 

Ignoring the EM fields, the elastic response of the neu- 
tral quantum hexatic can be straightforwardly calculated 
from Eqs. (0EJ) yielding similar results for the mode spec- 
trum as we discussed for the smectic Let us now 
focus on the electromagnetic response, the only way to 
probe the system from the outside given that we deal with 
electrons. According to standard linear response theory, 
the electrodynamic response is governed by the dielec- 
tric tensor, which is completely parametrized in terms 
of transversal and longitudinal dielectric functions. For 
isotropic media, these can be expressed in terms of diag- 
onal elements of the photon self-energy n as (2t| (with q 
the momentum and lu u the Matsubara frequencies), 



e L {q,uj n ) = 1 - nT ^" n) , e T {q,u n ) 



1 



(G) 



in the Coulomb gauge for the EM fields, d a A a — 0. We 
observe that the total dual action Eq. EJ) can schemat- 
ically be rewritten as, 
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where all the physics of the stress sector is lumped to- 
gether in the fully dressed stress gauge field propagator 
(G)ttvab- Given that the coupling between stress and EM 
gauge fields is linear, the photon self-energies are easily 
derived by straightforwardly integrating out the stress 
photons from Eq. J7J), 
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FIG. 1: The electromagnetic absorptions in the 'nearly or- 
dered' superconductor as seen by electron-energy loss spec- 
troscopy (Im[l/eL(q, to)]) as function of frequency ui and mo- 
mentum q. In this example, we take a typical Poisson ra- 
tion v = 0.28, characterizing the 'background' Wigner crystal, 
while the 'shear Higgs mass' fl — 0.05oj p and shear penetra- 
tion depth As = 1/qs = cr/fi are taken to be representative 
for the situation cuprate superconductors (uj p ~ 1 eV, Q ~ 50 
meV, As — 10 nm). Besides the strong plasmon pole dom- 
inating the long wavelength dielectric response, we also find 
a weak absorption which corresponds to the massive shear 
photon giving away the presence of the 'dual shear supercon- 
ductor', to be regarded as the unique fingerprint of a super- 
conductor characterized by transient translational order ex- 
tending over distances large compared to the lattice constant. 
In the inset, the weights of the plasmon (solid red, I p ) and 
the shear photon (dashed green, I s h) in the electron energy 
loss spectrum as function of momentum are given. The latter 
can be used to isolate the shear photon in the experimental 
data. At long wavelengths (q < qs) I 3 h grows quadratically 
with momentum, with a prefactor set by the ratio of the shear 
mass and the plasmon frequency 6 — Q/u) p , while it reaches 
a maximum at at intermediate wavelengths (qs < q < Qe)- 
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The coupling constants g a ,^v and the H bare are tabulated 
m Ref. UJ. It is a straightforward exercise to obtain 
explicit expressions for the dielectric functions Eq s. JSJl- 
As we already noticed in the earlier work |12| . from 
the transverse electromagnetic response it follows that 
the quantum hexatic is also an electromagnetic Meissner 
state: this should be surprising because this supercon- 
ductivity does not involve off-diagonal long range order 
of the constituent bosons and the presence of the dual 
dislocation state turns out to be a sufficient condition. 



In the remainder of this paper, we will focus on the lon- 
gitudinal electromagnetic response for which 
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(9) 



This is our main result: the longitudinal response is 
at least in principle accessible by experiment at the fi- 
nite momenta which are required to see the influences 
of the transient order: all that is needed is a measure- 
ment of the electron energy loss spectrum corresponding 
to Im[f/£i(g,w)] in the relevant kinematic regime. The 
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existence of two modes is evident from the poles. In the 
limit As/a — > oo, this expression has a universal status 
for small momenta; at larger momenta one has to account 
also for damping by perturbative mode couplings, single 
particle excitations, etcetera. In Fig. 1, we show an exam- 
ple of the loss-spectrum which could be representative for 
cuprates assuming that the characteristic disorder scales 
revealed by the measurements of the spin fluctuations are 
rooted in the charge fluctuations discussed here: as com- 
pared to the plasma frequency ~ 1 eV, the shear Higgs 
mass ft should be small, (~ 1CP 2 eV), while the charac- 
teristic length scale parametrized by As should be ~ 10 
nm), consistent with a typical electronic velocity ct — 1 
eVA. The bottom line is that besides the intense plas- 
mon one finds a weak feature at low energy in the loss 
spectrum, corresponding to the massive shear photon. 

The weight of the shear photon has a non-trivial de- 
pendence on both the ratio 8 = Q/u> p as well as on 
the momentum q. This is easily computed from Eq. 

and the results are summarized in inset of Fig. 1. 
In the 'liquid' long wavelength regime q < qs, the 
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weight of the shear photon I s h(o) — 4(1-32) 
I s h has a maximum at intermediate momenta q max — 
y/2(l — i/)/ (3(3 + u)) ~ 0.4g e where the shear photon ac- 
quires a weight relative to the plasmon {I s h{qmax) / Ip) — 
v) 8 2 which we expect to be of 
0.05: the elec- 



(3^3(1"- u))/(16y/3- 
order 0.1% in cuprates given that 9 
tromagnetic weight of the shear photon is expected to 
be very small. This momentum and 6 dependence of the 
electromagnetic strength of the shear photon can however 
be taken as its fingerprint, since it is rooted in the mech- 
anism through which the shear photon becomes electro- 
magnetically active. Both shear stress and dislocations 
forming the building blocks for the massive shear photon 
do not carry volume and are therefore electrically neu- 
tral. However, in the Wigner crystal and in the liquid, 
at distances small compared to As, the plasmon carries 
like longitudinal phonon shear components at finite mo- 
menta. These shear components have to be 'removed' 
from the plasmon, turning it into the purely compres- 
sional mode of the liquid at distances large compared 
to A s . This causes a linear mode coupling between the 
plasmon and the shear photon, with the latter 'steal- 
ing' some electromagnetic weight from the former. This 
explains the strong dependence on 6, the characteristic 
maximum, and the vanishing of the weight at q — > 
as even in the Wigner crystal, the plasmon turns into a 
purely compressional mode in the long wavelength limit. 

What has to be done to detect the shear photon? It 
is clear that one needs a very high resolution measure- 
ment of the loss function, given its low energy (~ 10 
meV) and very small spectral weight. Optical methods 
are not suited because these measure at momenta which 
can be regarded as infinitesimal even on the scale of qs- 
Obviously, the easy way is to measure directly the loss 



function in this kinematic regime. This requires a mil- 
lielectronvolt energy, and nanometer spatial resolution. 
Although such instruments are not available right now, it 
just seems to involve a relatively minor engineering effort 
to get into this regime: reflection (low energy) EELS has 
the resolution but for unknown reasons there are prob- 
lems measuring the electronic response |3l) ; transmission 
(high energy) EELS is limited to a resolution of order of 
~ 100 meV but there seems much room for improvement 
using modern electron optics. The strongest contender 
appears to be resonant inelastic soft X-Ray scattering, 
where at present a major instrumental development pro- 
gram is unfolding aimed at reaching meV resolution. 
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